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The Hubbard chain and spinless fermion chain are paradigms of strongly correlated systems, very well un¬ 
derstood using Bethe ansatz. Density Matrix Renormalization Group (DMRG) and field theory/renormalization 
group (RG) methods. They have been applied to one-dimensional materials and have provided important in¬ 
sights for understanding higher dimensional cases. Recently, a new interacting fermion model has been intro¬ 
duced, with possible applications to topological materials. It has a single Majorana fermion operator on each 
lattice site and interactions with the shortest possible range that involve 4 sites. We present a thorough analysis 
of the phase diagram of this model in one dimension using field theory/RG and DMRG methods. It includes a 
gapped supersymmetric region and a novel gapless phase with coexisting Luttinger liquid and Ising degrees of 
freedom. In addition to a first order transition, three critical points occur: tricritical Ising, Lifshitz and a novel 
generalization of the commensurate-incommensurate transition. We also survey various gapped phases of the 
system that arise when the translation symmetry is broken by dimerization and find both trivial and topological 
phases with 0, 1 and 2 Majorana zero modes bound to the edges of the chain with open boundary conditions. 

PACS numbers: 


I. INTRODUCTION 


Models of strongly correlated electrons have provided out¬ 
standing challenges to condensed matter theory for many 
decades. These models are expected to describe the physics 
of important materials such as high temperature supercon¬ 
ductors, complex oxides and quantum magnets. While ex¬ 
act or controlled analytical and numerical solutions are rare 
for the two and three dimensional cases, great understanding 
has been achieved for the one-dimensional (ID) versions of 
these models. This is true despite the fact that quantum fluc¬ 
tuations are generally enhanced in ID, leading to the break¬ 
down of Fermi liquid and mean-fleld theories. The remark¬ 
able progress may be ascribed to a number of theoretical tech¬ 
niques: the Bethe ansatz, numerical methods, notably the den¬ 
sity matrix renormalization group (DMRG) Q, and fleld the¬ 
ory methods based on bosonization and the renormalization 
group (RG). 

The canonical example of correlated electrons, namely, the 
Hubbard chain depicted in Fig. has interaction between 
electrons with opposite spin on the same site, 

H = + H.c.) + U (1.1) 

jo- j 

It has been well studied by all these methods and its phase 
diagram is understood in great detail for both signs of the on¬ 
site interaction U and for general values of the chemical po¬ 
tential p (which couples to the total density as -pTjjcr^j,cr)- 
(For a review see O.) In Eq. ( |I.I[ ) cjo- annihilates an elec¬ 
tron on site j with spin cr ='[,i and hjcr = jo-• An even 
simpler version of the Hubbard chain for spinless fermions, 
which we henceforth refer to as the Dirac chain, is described 


by the Hamiltonian 

// = X + H.c.) + V{hj - l/2)(n;+i - 1/2)1 • 

( 1 . 2 ) 

The model \\2) is illustrated in Fig.[2b). It has also been very 
well studied by all the methods mentioned above E m for 
general chemical potentials. In both models Coulomb interac¬ 
tions are treated as being highly screened, with strictly on-site 
interactions in the spinful Hubbard case and nearest-neighbor 
interactions, the shortest range possible due to the Pauli ex¬ 
clusion principle, for the spinless case (note that fp: = hj). 
The spinless fermion model is equivalent to the XXZ spin-1/2 
chain, providing important experimental realizations. 

There is a third very natural ID model which reduces the 
number of degrees of freedom per site by another factor of 
1/2 compared to the spinless fermion chain. The fermions in 
Eq. ^ are complex (Dirac) fermions with cj 4^ c\. The 
analog of Eq. for real (Majorana) fermion operators yj, 
obeying 

r] = r> {ry>ri) = (i-3) 

can be written as 

H = X + syjyj+\yj+2yj+3\ ■ (i -4) 

j 

The Hamiltonian again has nearest-neighbor hopping and in¬ 
teractions of the shortest possible range: Since = I, the 
shortest range nontrivial interaction spans four consecutive 
sites as shown in Fig. ^c). In this case, no chemical poten¬ 
tial term is possible so the model only has one dimensionless 
parameter g/t. Surprisingly, this third canonical model of in¬ 
teracting fermions has remained relatively unexplored EHH. 
This is partly due to the fact that Majorana fermions are yet 
to be observed as elementary particles of nature. However, 
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FIG. 1: (a) The Hubbard chain with on-site interactions, (b) The 
spinless Dirac chain with nearest-neighbor interactions, (c) The most 
local Majorana chain with 4-site interactions, (d) The phase dia¬ 
gram of the Majorana chain with the Hamiltonian ( |1.4| ) as a func¬ 
tion of g/t. Setting t = 1, it consists of four phases: a 4-fold de¬ 
generate gapped phase separated by a generalized commensurate- 
incommensurate (C-IC) transition at g ^ -2.86 from a critical phase 
with central charge c = 3/2 comprised of a critical Ising and a decou¬ 
pled Luttinger liquid (LL). The Ising-hLL phase is, in turn, separated 
from a critical Ising phase by a Lifshitz critical point with dynami¬ 
cal exponent z = 3 at g -0.28. For positive g, we have another 
transition from the Ising phase to a doubly degenerate supersymmet¬ 
ric gapped phase at g ^ 250 with the phase transition described by 
the tricritical Ising (TCI) conformal field theory (CFT) with central 
charge c = 7/10. 


Majorana zero modes (a special type of Majorana fermion oc¬ 
curring at exactly zero energy), are predicted to occur as col¬ 
lective degrees of freedom in several condensed matter sys¬ 
tems EHm, with remarkable experimental progress reported 
in recent years MM- 

The model defined by Eq. ( |1.4| ) above has been proposed to 
describe the low-energy physics of a ID vortex lattice formed 
in a superconducting film on the surface of a strong topolog¬ 
ical insulator (STI) 125^271 . A Majorana zero mode (MZM) 
exists in the core of each vortex rrn . The tunneling term 
in Eq. ( |1.4| ) makes these Majorana modes dispersive, mov¬ 
ing their energy away from zero (corresponding to an isolated 
MZM). By tuning the chemical potential of the STI, it is pos¬ 
sible to tune the hopping parameter t to zero (due to a chiral 
symmetry 125112811291), providing access to the strong cou¬ 
pling regime \g/t\ » 1. It was shown that ordinary Coulomb 
interactions in the STI lead to g < 0. However, the g > 0 
case might also be relevant due to the effectively attractive 
electron-electron interactions leading to superconductivity. 

Unlike the Hubbard or spinless Dirac fermion chain, the 
Majorana chain has no continuous symmetries and is not 
Bethe ansatz integrable, as far as we know. Also unlike the 


Hubbard and Dirac chains, the Majorana chain has a nontriv¬ 
ial ground state in the strong coupling limit, r = 0. To study 
the model numerically, and understand its behavior in certain 
limits, it is convenient to construct Dirac fermion operators 
from pairs of neighboring MZMs. Combining, e.g., sites 2 j 
with 27 -I- 1 into 


Cj = {y2j + iy2j+\)l2, (1.5) 

the Hamiltonian takes the form 

^ = Z f + O+i)] 

^ (1.6) 

+ g [-Pjpj+l + (c] - Cj)Pj+l(c]^2 + 0 + 2 )]} , 

where pj is shorthand for 2hj - 1. Note that the hop¬ 
ping term turns into a combination of a chemical poten¬ 
tial, nearest-neighbor hopping and nearest-neighbor pairing 
terms, whereas the interaction term turns into a combination 
of nearest-neighbor interactions and second-neighbor hopping 
and pairing with amplitude depending on the filling at the cen¬ 
tral site. As anticipated, negative g corresponds to a repulsive 
interaction term. Similar to the model of Eq. the Ma- 
jorana chain model can also be exactly mapped into a spin 
Hamiltonian by a Jordan-Wigner transformation cr^. = pj and 

cr+ _ which leads to 

j j j j 

(1.7) 

Despite the complicated form of the Hamiltonian ( |1.6| ) and 
the absence of particle number conservation, we have solved 
this model for large system sizes using DMRG as well as by 
a combination of field theory/RG and mean field arguments. 
We have determined the complete phase diagram for both 
signs of the interactions. Altogether we find 4 different sta¬ 
ble phases, sketched in Eig.[2d). Without loss of generality, 
we may assume ^ > 0. Eor -0.285 < g/t < 250, we find a 
critical “Ising” phase with a single gapless relativistic Majo¬ 
rana fermion excitation. Eor -2.86 < g/t < -0.285, we find 
a phase consisting of decoupled gapless Majorana and Lut¬ 
tinger liquid excitations, with central charge c = 3/2 and a 
continuously varying Luttinger parameter K < 1, correspond¬ 
ing to repulsive interactions. The Luttinger-liquid sector of 
this phase has a conserved charge, which is an emergent sym¬ 
metry, not present in the lattice model. At strong coupling of 
either sign, gapped phases with spontaneously broken sym¬ 
metries and degenerate ground states occur. Eor g/t < -2.86, 
the ground state is 4-fold degenerate. The symmetry of the 
Hamiltonian ( |1.4| ) under translation by one Majorana site is 
spontaneously broken down to translation by 4 Majorana sites 
in the ground-state wave functions, giving rise to a unit cell 
of 4 sites. Eor g/t > 250, on the other hand, the ground 
state is 2-fold degenerate with a unit cell of 2 sites. (Trans¬ 
lation symmetry is broken down to translation by two Majo¬ 
rana sites). Since two Majoranas can be combined into one 
Dirac fermion as in Eq. ( |1.6| ), in terms of Dirac fermions, 
translation invariance is not broken in the gapped phase with 

































3 


positive g (git > 250). However, the Hamiltonian ( |1.6[ ) is 
symmetric under the particle-hole transformation cj c\ at 
r = 0. This additional particle-hole symmetry, only present for 
^ = 0, is spontaneously broken in this limit. The models with 
gjt = +00 are equivalent and we find that for large positive 
git there is a low lying doublet of excited states, with energy 
per unit length cx \t\, which become degenerate with the two 
ground states as g/r ^ oo. 

There are 3 critical points in the phase diagram, as well as 
the first order transition dX git = oo. The transition from the 
Ising phase to the gapped 2-fold degenerate phase at git = 250 
is in the tricritical Ising universality class, with c = 7/10. The 
corresponding CFT is supersymmetric. The relevant operator 
that drives the transition at this critical point respects the su¬ 
persymmetry, which is therefore present in the gapped phase 
for git slightly bigger than the critical value git = 250. (The 
effective Hamiltonian is also supersymmetric at git slightly 
smaller than this critical value but supersymmetry is sponta¬ 
neously broken in this Ising phase.) On the negative-g side, 
the transition between Ising and Ising plus Luttinger liquid 
(LL) phases at g/r = -0.285 is characterized by the vanishing 
of the velocities, giving rise to a dynamical exponent z = 3 
at the critical point, and an effective Fermi wave vector which 
is zero at the critical point and grows as we move into the 
Ising-rLL phase. This critical point describes a Lifshitz tran¬ 
sition as the number of low-energy Majorana modes changes 
from one (in the Ising phase) to three (in the Ising-rLL phase), 
which is analogous to the change in the topology of the Fermi 
surface. 

The dynamically induced Fermi wave vector continuously 
changes within the Ising-rLL phase and reaches a commensu¬ 
rate value of 7r/4 at the third critical point, git = -2.86, where 
a transition occurs into a gapped 4-fold degenerate state. This 
is related to the commensurate-incommensurate (C-IC) tran¬ 
sition, which occurs in the spinless Dirac chain as the chem¬ 
ical potential is varied for V > 2t |[30l|3Tl. It is, however, a 
nontrivial generalization of the C-IC transition since the inter¬ 
action driving the transition couples LL and Ising sectors. We 
therefore term this transition generalized C-IC. 

Although the bulk of this paper deals with the transla- 
tionally invariant (translation by one Majorana site) system 
described by the Hamiltonian ( |L4| ), a more general one¬ 
dimensional array of interacting Majoranas, in which we only 
have invariance under translation by two Majorana sites, may 
be relevant to experiments. This may occur, e.g., if the vor¬ 
tices are arranged in a dimerized pattern. In general, the phase 
diagram of the dimerized model will depend on three dimen¬ 
sionless parameters gilti, and g 2 lti instead of just one 
parameter git [See Eq. (2.5)]. While an in-depth analysis 
of the phase diagram of the dimerized model is beyond the 
scope of this paper, we present a topological classification of 
the gapped phases of this model and identify them using self- 
consistent mean-field calculation. 

We argue that according to the results of Refs. |[32j|33l the 
gapped phases of the dimerized model fall into eight cate¬ 
gories characterized by the difference between the number of 
time-reversal even and odd Majorana bound states at the end¬ 
points of the chain with open boundary conditions. Within the 


self-consistent mean-field scheme, we find that four of these 
topological classes can (and do in parts of the phase diagram) 
appear. The mean-field picture also provides a good descrip¬ 
tion of certain features of the phase diagram of the nondimer- 
ized model. In particular, it accurately predicts the Lifshitz 
transition. 

The remainder of this paper is organized as follows. In the 
next section we discuss the weak and strong coupling limits 
(for both signs of the coupling constant, g). In Sec. Ill we 
analyze the model with g < 0 and in Sec. IV, g > 0. In 
Sec. V we employ mean-field theory to discuss a more gen¬ 
eral dimerized model, where t and g alternate, and present a 
topological classification of the gapped phases arising in the 
system.. Sec. VI elaborates on possible experimental realiza¬ 
tions of our model and various signatures of its phases and 
phase transitions. Conclusions are given in Sec. VII and fur¬ 
ther technical details appear in three Appendices. 

II. WEAK AND STRONG COUPLING LIMITS 
A. Weak Coupling Limit 

Weak coupling treatment of interacting Majoranas relies on 
solving the noninteracting problem and treating the interac¬ 
tions with RG (see Refs. |[34l - [36ll for a few examples).We first 
discuss the noninteracting case, g = 0. As shown in Appendix 
[A| the low-energy Hamiltonian of the noninteracting system 
is given by 

Ho = iv J" dx[yLdx'YL-yRdx7R], {v = At) (2.1) 

[jRldx), yRiLiy)} = (1 /2)(5(x - j), (2.2) 

which describes free relativistic Majoranas. This is a CFT 
with central charge c = 1/2, corresponding to the critical point 
of the transverse-field Ising model [see Eq. with g = 0]. 
Note that no mass term = im f dxyRyi appears. This is a 
consequence of the translational symmetry, yj yj+i which 
maps 

7r 7r, 7l -7l- (2.3) 

Taking into account that 7ril(^) vary slowly on the lattice 
scale, it is now straightforward to write the interaction term 
projected onto the low-energy states as 

Hint ^ -256g J" dxyRdxyRjLdxyL- (2.4) 

Note that 2 derivatives are necessary to get a nontrivial in¬ 
teraction by Eermi statistics. This interaction has RG scal¬ 
ing dimension 4, 1/2 for each Majorana field and 1 for each 
derivative, and so is highly irrelevant. [Dimension A > 2 is 
irrelevant in relativistic (1 + 1) dimensional field theory.] 

This implies the existence of an extended phase with a 
massless Majorana in the vicinity of g = 0. Note that the 
present scenario is quite different from what happens in the ID 
Hubbard and Dirac [see Eq. ( |L2| )] chains, where nonderiva¬ 
tive operators can occur. These lead to the Mott transition at 
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infinitesimal coupling in the Hubbard chain and to the contin¬ 
uously varying Luttinger parameter in the gapless phase of the 
spinless Dirac chain. 

B. Strong Coupling Limit 

Unlike the spinless Dirac chain, the strong coupling ground 
states of the Majorana chain are nontrivial. However, in¬ 
sight into the nature of these ground states can be obtained by 
considering a model with alternating hopping and interaction 
terms: 

^ “ X + ihyij+iyij+i 

j 

+ ^1 y 2 jy 2 j + iy 2 j + 2 y 2 j +3 + ^2727+1727+2727+3727+4) • 

(2.5) 

We will argue that the symmetry of translation by one site is 
spontaneously broken at strong coupling and that the ground 
states are thus qualitatively similar to the ones for the Hamil¬ 
tonian of Eq. 

In the strong coupling limit ti = 0, the model is trivially 
soluble if either gi or g 2 = 0. For instance, when g 2 = 0 it 
follows from Eq. ( |1.6| ) that 

H ^-giYjPjPj+i- ( 2 . 6 ) 

j 

For > 0, corresponding to attractive interactions, the two 
ground states have all Dirac levels filled (pj = 1) or empty 
(Pj = -1). The case of gi >0 corresponds to a ferromagnetic 
Ising interactions in the spin representation [see Eq. O]- 
For gi < 0, a charge density wave occurs, with every second 
Dirac level filled or empty, pj = ±(-l)^ (antiferromagnetic 
spin chain). Similarly, if gi = 0, we combine Majoranas to 
form Dirac operators on sites {2j + 1,27 -r 2) obtaining the 
same Eq. \2.6) for a different set of Dirac fermions (shifted 
by one Majorana site). These 4 different ground states are 
indicated in Figs. and for the cases gi > 0 and gi < 0 
respectively. Clearly, there is a gap, 2|g/|, to the lowest energy 
excited states. While we can only solve the model exactly 
at strong coupling when the interactions vanish on half the 
quartets of Majoranas, we will see that the 4 corresponding 
ground states, for either sign of gi, are cartoon representations 
of the actual ground states of the uniform chain in the strong 
coupling limit. If this is correct, there are two spontaneously 
broken discrete symmetries. One of these is translation by 1 
site: jj yj+\- In the special case t = Q, the model has 
particle hole symmetry, cj c \. For both definitions of Dirac 
modes this corresponds, up to a phase, to 

77 ^ (-1V77- (2.7) 

Clearly this particle-hole symmetry is spontaneously broken 
when ti = t 2 = ^ and either gi or g 2 = 0. The two ground 
states in those cases are mapped into each other by a particle- 
hole transformation. We expect that it is also spontaneously 
broken when gi = g 2 and t\ = t 2 = 0. The combination of 


spontaneously broken translation symmetry and particle-hole 
symmetry results in 4-fold ground state degeneracy. 

It is also important to note that, for r = 0, the duality trans¬ 
formation: 

yj -Jp (; = 4 «) 

^ yj, (otherwise) (2.8) 

maps the Hamiltonian H -H corresponding to g ^ -g. 
This transformation interchanges the cartoon ground states of 
Figs. © and 0. Thus if the spontaneously broken ground 
states occur for one sign of g they must occur for the other. 

We note that the Hamiltonian ( |2.5| ) can be exactly mapped 
into a spin Hamiltonian by a Jordan-Wigner transformation, 
giving 

//=F - ^2 F Z - ^2 z 

j j j j 

(2.9) 

The special case gi =0 corresponds to the anisotropic 
next-nearest-neighbor Ising (ANNNI) model. (For a re¬ 
view see E3-) For ti = t 2 = 0, the Hamiltonian 
maps to a multispin Ising chain with H = -giYuj^) ~ 
g 2 Yuj^]-\^]^)+i ^)+2 through a transformation r^. = 

This model has been studied in the literature |[38ti4ll and 
based on a connection with the 8-state Potts model as well as 
some numerical evidence iOlED is believed to have a first- 
order transition at the self-dual point gi = g 2 . Therefore, (i) 
the system is likely gapped for ti = ^2 = 0, gi = g 2 and (ii) the 
ground states on two sides of the transition (respectively dom¬ 
inated by gi and g 2 only) are present at the self-dual point. In 
the proceeding sections, we provide numerical evidence for 
this scenario. 

Assuming these gapped broken symmetry ground states at 
^ = 0, let us now consider the likely effect of a small nonzero 
t. It is again convenient to start with the case g 2 = 0, gi = g 
and turn on only a small ti in the dimerized Hamiltonian of 
Eq. (|g. Again, in Dirac notation, we have a special case of 
Hamiltonian ( |1.6| ): 

^ = Z (^‘^2 - SiPjPj+\) ■ (2.10) 

j 

We see that t breaks particle-hole symmetry, favoring empty 
Dirac levels for t\ > 0 or filled for i <0. Its effects are quite 
different depending on the sign of gi. For gi > 0, ti splits the 
degeneracy between the two ground states, favoring the state 
with all Dirac levels empty or filled, depending on its sign. 
This would correspond to a first-order phase transition with a 
jump in (pj) (and consequently a jump in the occupation num¬ 
ber (hj) as Pj = 2hj - 1) at =0 (see Appendix C for nu¬ 
merical evidence). On the other hand, for gi < 0, turning on a 
small ti does not split the degeneracy of the two ground states 
which remain the exact ground states for l^il < |gi|, at which 
point a transition occurs (through a multicritical point (631) 
into one of the gi >0 ground states with all Dirac levels empty 
or filled, depending on the sign of ti . The Ising spin-chain rep¬ 
resentation Pj = cr^. is also illuminating: the term proportional 
to ti in Eq. ( |2.10| ) serves as a magnetic field coupled to the 
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^1 > 0 



- # 

^2 > 0 

- # 



FIG. 2: Mean-field ground states for g > 0, corresponding to ferro¬ 
magnetic states in the spin representation. The small blue and green 
circles represent the sites of the Majorana chain. Bold links with 
large circles on them represent Dirac fermions formed by combin¬ 
ing two Majoranas (two types of combinations are considered: blue- 
green and green-blue). A filled circle corresponds to the Dirac level 
being filled and an empty circle to it being empty. 


magnetization, which splits (does not split) the degeneracy of 
a ferromagnet (antiferromagnet) when it is less than the ex¬ 
change coupling. Note that, in this sub-section, we have ex¬ 
plicitly broken the symmetry of translation by one Majorana 
lattice site. Given this assumption, we have argued that, for 
g < 0, but not g > 0, there is a further spontaneous breaking 
of symmetry of translation by two Majorana sites. We will 
argue that these results are qualitatively correct for the trans- 
lationally invariant Hamiltonian, with t\ = t 2 , gi = g 2 , In the 
strong coupling regime, \g\ » \t\, there is a 2-fold degener¬ 
ate ground state for large positive g and a 4-fold degenerate 
ground state for large negative g. However, we shall see that 
the transition out of the strong coupling phases is continuous, 
for both signs of g. 


III. REPULSIVE INTERACTIONS: g<0 


As discussed in Sec. EE we expect the Ising phase to per¬ 
sist up to a finite critical g/t of either sign. We have con¬ 
firmed this primarily by calculating several energy levels in 
the finite-size spectrum, with antiperiodic boundary condi¬ 
tions (APBC). We can classify all states by their fermion par¬ 
ity which we formally define, for a chain with an even number 


< 0 




♦ 


--•-0- 

^2 < 0 

-•-0-- 

FIG. 3: Mean-field ground states for g < 0, corresponding to the 
antiferromagnetic states in the spin representation 



FIG. 4: Velocity in the Ising phase near the Lifshitz transition. As 
the transition is driven by a renormalization of the dispersion relation, 
mean-field calculations are in approximate agreement with DMRG. 


of Majorana sites L = 2l numbered 0,1,2,... (L - 1), as 


e-i 

F = Y[{}y2iy2M)- ( 3 - 1 ) 

Note that, defining Dirac fermions by Eq. ( |1.5| ), this becomes 
F = nfjo Pj - where Nf is the number of Dirac 

fermions added to the vacuum state. [While Np is not con¬ 
served, (-1)^^ is in Hamiltonian ( |1.6| ).] In the Ising, i.e., 
free Majorana, phase, the low-energy excitations with APBC 
simply correspond to creating particles of energy £’v|^|, with 
k = 2nv{n -f l/2)/L for integer n. Thus the lowest energy 
excited state has opposite fermion parity to the ground state, 
energy ttv/L and is two-fold degenerate. The lowest excited 
state with the same fermion parity as the ground state has 2 
fermions added at ^ = tt/L and k = n - nlL with energy 
2nvlL (k is restricted to the interval 0 < ^ < tt. See Ap¬ 
pendix A.). This ratio of excitation energies, two, remains 
essentially constant in our DMRG results over the entire Ising 
phase while the velocity, v varies. We performed the computa¬ 
tions at g = 1 by varying t. As tig -3.512 (corresponding 
to gjt = -0.285), we find that the velocity vanishes as shown 
in Fig. 1^ As we decrease t, the behavior of the finite-size 
spectrum changes (compared with the Ising phase realized at 
larger tjg). As we see this is a Lifshitz transition correspond¬ 
ing to a change in the topology of the Fermi surface. The ve¬ 
locities were extracted from a linear fit of the finite-size gaps 
computed with DMRG to 1 /L for a range of system sizes with 
L = 40,... 200 Majoranas. We kept 500 states in the DMRG 
computations to achieve convergence in the energy gaps. 

Computing the central charge sheds light on the nature of 
the phase for larger hopping t. In a critical phase described by 
a (1 -f 1)J CFT, the entanglement entropy 5 of a subsystem 
of length y with the rest of the system is related to the central 
charge c through 


5 = 



-f const. 


(3.2) 


for a system of length L with periodic boundary condi¬ 
tions m. We numerically computed the entanglement en¬ 
tropy with DMRG (keeping 300 states) and extracted the cen¬ 
tral charge from the relationship above. As expected, the cen¬ 
tral charge in the Ising phase is c = 1/2. However, as seen in 
Fig. it jumps to c = 3/2 at the same value of t, for which 
the velocity goes to zero. This indicates that upon decreas¬ 
ing the value of t beyond this transition point, three species of 
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FIG. 5: The behavior of the central charge as a function of t for 
g = -1- 


low-energy Majoranas appear at this phase transition. In fact, 
as discussed below, the same behavior can arise in a noninter¬ 
acting model with third-neighbor hopping. Therefore, we can 
understand this phase transition simply in terms of a kinetic 
energy renormalized by interactions. Incidentally, as we will 
see in Sec.jVj a mean-field calculation captures this transition 
with good accuracy (see Fig.[^. 

Third-neighbor hopping is indeed allowed by all symme¬ 
tries. [Spatial parity symmetry yj (-lyy-j forbids a sec¬ 
ond neighbor hopping term. Notice that a naive parity trans¬ 
formation yj y-j changes the sign of the nearest-neighbor 
hopping term and the (-1)^ term simply correct for this.] Con¬ 
sider a quadratic Hamiltonian 

H = iY, + f'ri+3] = :^ y Eky{-k)y{k) (3.3) 

j ^ k 

with Ek = -f Af sin(3^). As in Appendix it is con¬ 

venient to regard y{k) as an annihilation operator for the re¬ 
gions of k where > {) and write y{k) as for the 

complementary regions. Consider the case ^ > 0, < 0. 

For f > -t/3, Ek vanishes at ^ = 0 and n only, with ve¬ 
locity V = At + \2t'. However, v ^ 0 at = -^/3. For 
t' < -tl3, Ek vanishes at 4 other points, ±ko and +(;: - ^o) 
with sin^o = (1/2) ^/3 -f tjf. Now there are 3 regions of k 
where E^ > 0, shown by thick black lines in Fig. ([^. The 
velocity at ^ = 0 is vq = Ibsin^^o, while at ^ = ko^n - ko, 
we have v = 2vocos^o- Note that vq and v increase linearly 
with —t' — tl3 while k^ increase more rapidly cx V-f - f/3. 
Here, k^ plays the role of a Fermi wave vector. We may again 
introduce relativistic fermions to represent the low-energy ex¬ 
citations as 

yj ^ 2rL(7)+(-l)-'27K(7)+[e-'^»V«(;) + Mj) + 

Here iI/r/l are Dirac fermion operators, simply related to the 
Fourier modes of the original Majoranas as 

= 7(ko + q), ^d-q) = y(7i: - ko - q), -A < q < A, 

(3.5) 

where A 1 is the momentum cut-off of the low-energy sec¬ 
tor. Note that the right/left movers occur at k points where 
Ek has positive/negative slope. For k slightly larger than ko, 
y(k) is identified with a right-moving particle annihilation op¬ 
erator whereas for k slightly less than ko it is identified with 



FIG. 6: Dispersion relation indicating a Lifshitz transition. For \t'\ < 
t/3, the zeros of the dispersion relation are at ^ = 0 and k = tt. At 
\t'\ = t/3 the velocities of both these low-energy modes vanish and 
for \t'\ > t/3, new low-energy modes appear at finite momenta ko and 
n - ko, with ko - 0 dii the Lifshitz transition \t'\ - 113. 


a right-moving antiparticle creation operator. The low-energy 
Hamiltonian becomes 

//o = / I dxivQiyid^yL - yRd^yR) + 

(3.6) 

We now consider the effect of the interactions. These are 
most rigorously treated if we added a f term to the Hamil¬ 
tonian by hand, and then turned on a small g. However, we 
expect the universal properties of the resulting phase to also 
describe the case at hand where f is generated dynamically. 
But, in this case we are not in the weak coupling regime since 
\g\ must be 0(t) to drive the Lifshitz transition. Since we have 
more fields in the low-energy field theory, it is possible to have 
nonderivative interaction terms. Many of these come with spa¬ 
tially oscillating factors, making them irrelevant for general 
values of k^. However, there are two nonoscillatory 4-fermion 
interactions allowed by symmetry: 

Hint ~ I dx [go : ■ +g'yRyL (lAz-lAs + > 

(3.7) 

where go = -16g[cos “Cos(3^o)] for weak coupling and 
indicates normal ordering. Since we are considering g < 0, 
we have go > 0 corresponding to repulsive interactions. The 
effects of this term by itself are well-known and easily treated 
using bosonization techniques, leading to a Luttinger liquid 
(LL). This corresponds to a free massless relativistic boson 
theory with the RG scaling dimensions varying continuously. 
These scaling dimensions are controlled by a single dimen¬ 
sionless parameter K known as the Luttinger parameter, which 
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takes the value K = 1 - ^ for weak coupling. Generally, 
we have ^ < 1 for repulsive interactions. 

We now argue that the second term in Eq. ( |3.7| ) above is 
irrelevant in the RG sense for ^ < 1. The scaling dimension 
of il/L^R appearing in the g' interaction is IjK, which is larger 
than one for repulsive interactions. The yRjL factor in this 
term also contributes 1 to the scaling dimension which leads 
to A = 1 + 1/^ > 2, making it irrelevant. Therefore, we 
find that the Ising and LL sectors are decoupled in the low- 
energy theory. This implies that a f/(l) charge conservation 
symmetry emerges in the LL sector of the Ising-hLL phase, 
along with an associated Lermi wave-vector k^. The charge 
is related to the occupation of modes near momenta +ko and 
±{7i - ko) as 

N= L + q)y(h + q) + yHn - h - q)y(n - ka - q) ■ 

-A<q<A 

(3.8) 


by 1/ ^/K and 6 by ^/K. It then follows that Eq. (3.10) be 


comes 


AE 


In 

T 


AK 


(N - koL/ny 


vK 
+ — 
4 




(3.11) 

There is no change in the energy of Ising excitations due to 
the interactions, since they only act in the LL sector. 

The finite-size spectrum now exhibits great complexity as 
koLjn is varied, either by varying L or ko via the hopping pa¬ 
rameter t (with g held fixed at g = -1). We are interested in 
the ground state and the first excited state in the two sectors 
corresponding to the total fermion parity. While the absolute 
ground state must have Nj = Mj = 0, we may have Ising exci¬ 
tations in the ground state in a given parity sector. To proceed, 
e first define /even/odd (•^) as the closest even/odd integer to a 
eal number x. It is easy to show that these functions are given 
by 


■E 


The constant -1 in the sum above is chosen such that there are 
no fermions in the ground state. Lor example, instead of writ¬ 
ing the number of right movers as 2 ]_a<^<a which 

is nonzero in the ground state (with q < Q modes occupied), 
we write 

0<q<A -A<^<0 


/even(x) = ^ (2LxJ + 1 - (-1)^) , (3.12) 

/odd(x) = I(2LxJ + l + (-l)W), (3.13) 

where [vj is the floor function of v. 

Now, fixing the fermion parity, we can write the candidates 
for the ground state energy in each sector as 


which measures the charge of the right-movers with respect 
to the Lermi level. The above expression and its analog for 
the left-movers leads to Eq. ( |3.8| ) through N = Nr Nr [see 
Eq. (|^]. 

All of this is in good agreement with our DMRG results. 
Again, the finite-size spectrum provides a powerful technique 
for confirming the phase diagram and extracting the univer¬ 
sal parameters, which govern it. As shown in Appendix 
the excitation spectrum of the system for a noninteracting LL 
(^ = 1) can be computed using elementary methods and is 
given by 



^(A^ - koL/nf + 


, (3.10) 


where N and M (Nj and Mj) are integers of the same parity, 
labeling excitations in the LL (Ising) sector. Physically, N 
is the total number of particles (both right and left movers), 
while M is the difference between the number of right and left 
movers, characterizing the charge current (this also implies 
that N and M have the same parity). In the ground state, we 
either have M = 0 and N is the closest even integer to k^Ljn 
or M = +1 with N is the closest odd integer. The value of N 
in the ground state determines how much it is shifted from the 
eigenvalues of the operator N in Eq. p.8| ). We have neglected 
the particle-hole excitations, which do not contribute to the 
low-energy excited states we study in this paper. 

The effect of the marginal coupling constant go in Eq. ( [3.7| ) 
on the finite-size spectrum can be accounted for as follows. 
After bosonizing, the charge density becomes proportional to 
the spatial derivative of a massless boson field, dxf. Simi¬ 
larly, the current density becomes proportional to the deriva¬ 
tive of its dual boson, dxO. The Luttinger parameter scales f 


reven 
^0,Dirac 


rcven 
^0,Ising 


rodd 
^0,Dirac 


rodd 
^0,Ising 


2n 

2n 

2n 

Yi 

2n 


U.n(-L) - -L 
n n 


l2 


(3.14) 


n n 

kn kn 

fodd(-L)--L 

n n 

fe.en(-L) - -L 
n n 


l2 


nKv nvo 

-r--r —,(3.15) 

2L L ^ ^ 

(3.16, 


l2 


TTVo 


(3.17) 


where the superscripts indicate the total fermion parity and 
the subscript Dirac (Ising) denoted the absence (presence) of 
an Ising excitation in the ground state of the given sector. We 
then have 


E™ = nun(£-“,,,£™p, (3.18) 

(3.19) 

The absolute ground state energy is given by 
min(£g''“,£°id) = nun(£^;“^^, When is 

close to an even integer, the ground state has even fermion 
parity with N even and M = 0. However, if k^Ljn is close 
to an odd integer, the two-fold degenerate ground states have 
odd fermion parity with N an odd integer and M = +1. As 
K gets smaller, the latter scenario occurs for a larger range of 
koLIn determined from the condition 

< -(-1)L^«^/'^J {koLIn -\koLln\ - 1/2). (3.20) 

V 

The lowest excited state in the odd sector exhibits even 
greater complexity. Lor some ranges of parameters it can be in 
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FIG. 7: (a) Odd-even and even-even gap, scaled by L for g = -1 
and t - 2.25. The stars are DMRG data and the other symbols are 
the best fit to the LL-hlsing spectrum. The blue symbols correspond 
to excitations in the LL sector whereas the red ones correspond to 
excitations in the Ising sector. The black circle corresponds to an 
excited state of even parity containing simultaneous Ising and LL 
excitations, (b) The Majorana correlation function. The fit to spec¬ 
trum give K - 0.4517 and - 0.5444, while fitting the correlation 
function (for even x so that there are no (-1)^ oscillations) gives 
K - 0.4611 and Uq - 0.5375, in very good agreement. The DMRG 
results were checked for convergence (with negligible truncation er¬ 
ror) for each system size. The maximum number of states kept in the 
computations was 700. 


the Ising sector, with Nj = Mj = \ and for other ranges in the 
LL sector. This can be determined by comparing the energies 
of several candidates as in the ground-state case. A sample of 
such data, for t = 2.25 is shown in Fig.[^ We can fit the data 
using the 4 parameters vq, v, and K. As expected, the fit 
gets better as L increases. 

We also verified that the fermionic correlation functions are 
in agreement with the spectrum. We expect from Eq. p.4| ) that 

/(7o7jc+i) ~ + b sin (kox 0) (3.21) 

for even x, nonuniversal coefficients a and b, and phase 
shift 0. Using even x simplifies the fit as further oscil¬ 
lations with wave vector tt and tt - ko, i.e., (-ly/x and 
(-ly sm(kox (p') [see Eq. p.4| )], do not appear 

for even x. The values of K and ko extracted from this are in 
agreement with those obtained from the finite-size spectrum. 

We extracted the 4 fitting parameters from the finite-size 
spectrum and the results are plotted versus t in Eig. Note 
that both V and vq appear to increase linearly with ti-t while 
ko increases as - t as expected from the noninteracting 
model, ti is the value of t at the Lifshitz transition for g = 
-1. Also K decreases with increasing \g\lt as expected for 
increasing repulsive interaction strength. 



t 


FIG. 8: Parameters of LL-rising phase determined from fitting the 
even-even and even-odd gaps, for a range of t, with g = -1. 


As we decrease t to around 0.5, the quality of fits decreases. 
We expect a phase transition to a gapped phase for some tc < 
0.5. Importantly, it appears that ko 7r/4 and K 1/4 as 
we approach the phase transition and fitting to the Ising-fLL 
begins to fail. This transition can be understood by observing 
that a possible interaction term, which was not included in 
Eq. ( |3.7| ) due to its fast oscillations at generic ko, namely 

(X J dxyRjL - H.cJ , (3.22) 

becomes nonoscillatory at ko = njA. The scaling dimension 
of H' is I 4K. Our DMRG data in Eig. indicates that 
K 1 /4 at the transition point so this interaction becomes 
relevant. Since H' couples together LL-rising fermions, we 
expect to gap out both sectors. The fact that its dimension 
goes to 2 at precisely the point where it becomes commensu¬ 
rate appears to be a novel generalization of the commensurate- 
incommensurate transition, which occurs in the spinless Dirac 
chain MM- 

To determine the critical value of tc for the C-IC transi¬ 
tion, we utilize an exact degeneracy of the gapped phase with 
periodic boundary conditions (PBC) as discussed below. Eor 


applies, with a 4-fold degenerate gapped ground state. The 
gap appears to remain very small all the way to r = 0 mak¬ 
ing this phase challenging to study with DMRG. However, 
we found that the 4-fold ground state degeneracy provides a 
convenient way of determining tc accurately. This degener¬ 
acy turns out to be exact for L = 8A sites with PBC, where 
all states come in degenerate pairs of opposite fermion parity. 
This follows, for periodic boundary conditions (PBC), since 
translation by one site, T, maps F into -F: 

F = Ton • • • r 2 L-i ^ 7172 - - - 72L-170 = -F. (3.23) 

It then follows that for any energy eigenstate |0), T\i//} is also 
an eigenstate of the same energy and opposite fermion parity. 
The full 4-fold degeneracy of the ground state is less generic 
and signals the broken symmetry phase. It corresponds to 2- 
fold degeneracy of the ground states in both even and odd 


smaller t, we expect that the qualitative description in Sec. IIB 
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fermion parity sectors. Let us focus, for example, on the even 
fermion parity sector. It is then convenient to start with the 
case ^2 = ^ = 0, so that the Hamiltonian only contains the gi 
term. Then the two ground states are known exactly: 

= c|c^ ...4^_j|0> 

|^^ 2 > = 44 ... 4^10). (3.24) 

Now consider translation by 2 sites (in the Majorana chain), a 
symmetry even when only 9 ^ 0 : 

= Wi) 

= 44 ■ • • ( 3 . 25 ) 

Thus we can construct linear combinations of these degener¬ 
ate ground states: 

|-A±) = (|^^i)t#2))/3^. (3.26) 

which are eigenstates of with eigenvalues +/. This cor¬ 
responds to momentum ±71 j2 with respect to translation by 
2 Majorana sites, corresponding to 1 Dirac site. That is, we 
define momentum P' in this case by . Now con¬ 

sider gradually increasing g 2 to g\. Although the 2 ground 
states become more complicated they must remain degenerate 
since they have equal and opposite momentum, by spatial par¬ 
ity symmetry. At g 2 = ^1 the degeneracy becomes 4-fold due 
to the additional symmetry of translation by one site. Now 
defining momentum by T = , we can form linear super¬ 

positions of even and odd fermion parity ground states with 
momentum + 7 r/ 4 , + 37 r/ 4 . Turning on a small t, the ground 
state degeneracy must survive up to the phase transition at tc 
by the same argument. On the other hand, if the number of 
Majorana sites is L = 8 A -h 4 this argument fails, and mixing 
and splitting of the 2 low-lying states in each fermion parity 
sector occurs for finite system size. 

We then determine tc (for g = -1 ) from our DMRG data 
by finding the largest value of t at which there is an exact 4- 
fold ground state degeneracy (exact 2 -fold degeneracy with, 
say, even fermion parity) for L = SN and extrapolating this 
value of t to large £ = Ljl. The results are shown in Fig. 
We plot the value of tc, where this exact degeneracy splits as a 
function of 1 jt. Linear {a/i+tc) and quadratic (a/^-\-h/£^-\-tc) 
fits to 1 /£ give the same result for the extrapolation within the 
error bar of 0.025 and we thus find tdg = -0.35 + 0.025 
{gitc = - 2 . 86 ). 


IV. ATTRACTIVE INTERACTIONS AND THE 
TRICRITICAL ISING POINT: g > 0 

The system for g > 0 was studied in Ref. |43l where it was 
shown to realize the tricritical Ising model, a supersymmetric 
CFT with central charge c = 7/10, at a critical value oi git, 
which separates the c = 1/2 Ising phase (at weak coupling) 
from a doubly-degenerate gapped phase (at strong coupling). 
While the primary focus of this paper is on repulsive interac¬ 
tions, in this section, we review some of the salient results of 
Ref. [43] for completeness. We also analyze a first-order phase 




1 /^ 


FIG. 9: (a) The gap from the ground state to the first ecited state 
as a function of the hopping amplitude t for a fixed system size 
and fermion parity. To distinguish exact degeneracy (in the gapped 
phase) from a small gap (in the critical phase in the vicinity of tc), 
highly accurate DMRG computations were performed by retaining 
up to 1000 states. (b)The extrapolation of tc for the C-IC transition 
gives tc ~ 0.35, setting g = -1. 



FIG. 10: The gap ratio as a function of system size computed with 
DMRG (keeping 600 states) . A plateau emerges at tjg ^ 0.004 at 
the expected value of 7/2 for the TCI CFT. 


transition that occurs between the symmetry broken phases at 
g/t = 00 (fixed g > 0 and t = 0 ), which has not been discussed 
elsewhere. 

The arguments of Sec. |II A| supporting a critical Ising phase 
around g = 0, were independent of the sign of g. Therefore, 
the c = 1 /2 Ising phase is expected to extend to a finite value 
of interaction strength g also for attractive interactions. As we 
argued in Sec. B unlike the 4-fold degeneracy of the strong¬ 
coupling limit gapped phase for g < 0 , the gapped phase at 
strong coupling is doubly degenerate in the case of g > 0. A 
priori, it is not obvious that there is only one phase transition 
between these two strong- and weak-coupling phases. If this 
is the case, however, it is well known from the theory of an 
Ising model with vacancies 1441 . that the most natural phase 
transition between the critical Ising phase and a doubly de¬ 
generate gapped phase is the TCI CFT. This CFT is of great 
interest as it provides a rare example of emergent supersym¬ 
metry in condensed matter physics. 

It turns out that the doubly degenerate gapped phase has a 
very large correlation length. Therefore, numerical verifica- 
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tion of the above-mentioned scenario is exceedingly difficult. 
As showed in Ref.|43j however, universal ratios in energy gaps 
provide a powerful method of verifying the scenario and iden¬ 
tifying the transition point corresponding to the TCI CFT. Our 
primary diagnostic is the ratio 


R = 


zT-odd _ iT'even 

IT'even _ reven ’ 
^1 ^0 


(4.1) 


The relevant Hamiltonian can be conveniently rewrit¬ 
ten in terms of new Majorana operators aj = 727 and JSj = 
727+1 2 iS H = Hq H' with 

Ho = i {haj/3j + t^PjOj+i) (5.1) 

j 

H ^ g2^7^7+1^7+1^7+2^ • (5.2) 

j 


where the superscripts even and odd indicate the fermion par¬ 
ity, Eq (El) is the ground (first excited) state energy in the 
parity sector, and all energies are computed with antiperiodic 
boundary conditions. The finite-size spectra of CFTs can be 
derived from the scaling dimensions of their primary opera¬ 
tors. In particular, for the TCI CFT we have 


^Tci = 7/2. 


(4.2) 


As shown in Fig. 10 the numerically computed value of R 
plateaus at tig = 0.004 as a function of system size precisely 
at /? = 7/2 predicted for the TCI CFT. 

Several other universal ratios support the presence of a tran¬ 
sition between the Ising phase and the gapped phase through 
the TCI CFT. Unlike the c = 1 /2 Ising model, the low-energy 
fermionic excitations of the TCI CFT are not regular free Ma¬ 
jorana fermions. This indicates nontrivial exponents for the 
Green’s function {y(x, t)y(0, 0)). In particular, the equal-time 
Green’s function decays as ^ for the TCI CFT as opposed to 
- for the Ising case, as was verified in Ref. [43l 


V. SELF-CONSISTENT MEAN FIELD THEORY AND 
TOPOLOGICAL CLASSIFICATION OF THE GAPPED 
PHASES 


In this section, we extend the mean-field-like picture of 
Sec. |IIB| which provided a qualitative description of the 
strong-coupling limit, to a more systematic self-consistent 
mean-filed treatment. As argued in Sec. Ill the Lifshitz tran¬ 
sition can be understood in terms of the generation of an ef¬ 
fective third-neighbor hopping in the renormalization group, 
which can renormalize the dispersion relation. Such terms 
are also naturally generated in the mean-filed decomposition, 
and as we show, the self-consistent mean field theory can pro¬ 
vide an accurate description of this transition. We also pro¬ 
vide a topological classification of the more general dimerized 
model defined in Eq. ( |2.5| ). 

Many ID systems, such as poly acetylene, are known to 
spontaneously dimerize. In our present model, as realized by 
the ID chain of vortices, dimerization may also occur, lead¬ 
ing to a model with explicitly broken translational symmetry. 
Whether or not such a dimerization occurs and its amplitude 
will depend on the details of the vortex lattice physics, specif¬ 
ically the intervortex interactions and vortex pinning. We do 
not attempt to specify the conditions under which dimeriza¬ 
tion may take place or calculate its strength. Rather, we ex¬ 
plore in this Section the phases of the dimerized model and 
discuss their topological properties. 


We continue assuming that ti and t 2 are positive but con¬ 
sider either sign of gi and g 2 - Although the original trans¬ 
lation symmetry by one Majorana site is broken the model 
obeys the antiunitary time-reversal symmetry T generated by 
((Xj,l3j) (a j, -jSj) and i -i. We note that in the lan¬ 

guage of previous Sections this corresponds to yR ^ y^ 
and i which is the proper form of time-reversal in 

the relativistic field theory. In the absence of interactions it 
also obeys the antiunitary particle-hole duality C generated by 
((Xj,l3j) (^j^Pj) and i -i which maps Hq -Hq. This 

puts the noninteracting Hamiltonian Hq into the BDI class un¬ 
der the Altland-Zirnbauer classification EH. In ID its gapped 
topological phases are therefore classified by an integer in¬ 
variant y (461. If we define Va/js as the number of unpaired 
MZMs of type a IP bound to the left end of the chain with 
open boundary conditions then the invariant v coincides with 
Va - yj 3 . Note that the two types of Majorana modes are dis¬ 
tinguishable because they transform as even and odd, respec¬ 
tively, under 7“. Also, because there is the same number of a's 
and P's in the system for each a bound to the left edge there 
must be one p bound to the right edge and vice versa. 

In their seminal work Fidkowski and Kitaev (321 [33l 
showed that when interactions preserving T are added to such 
a Hamiltonian, which is the case here, its integer classification 
is changed to the Zg classification: phases characterized by in¬ 
variants y and y -r 8 become indistinguishable. 

It is important to note that symmetry T of Hq hinges on 
the absence of the second neighbor tunneling because terms 
such as iajaj+i would clearly break T. In a physical system 
this may be realized to a good approximation due to the ex¬ 
ponential decay of the MZM wavefunctions. Alternately, as 
demonstrated in Ref. (26l[27l, the T symmetry can be imple¬ 
mented exactly when the chain is realized in a system of alter¬ 
nating vortices and antivortices in the surface of an STI with 
the chemical potential tuned to the neutrality point. When T 
is broken the noninteracting system is in symmetry class D 
and its classification in ID is Z 2 , with or without interactions. 
Physically, then, only phases with even and odd index y are 
distinct. 

The Hamiltonian H contains 3 dimensionless parameters. 
A detailed analysis of this 3-dimensional parameter space us¬ 
ing DMRG and exact diagonalization would require a lot of 
computer time and we leave it to future studies. Here, we per¬ 
form a survey of its gapped phases using the mean-field the¬ 
ory. Although such MF theories often fail to accurately cap¬ 
ture the physics of systems in low dimensions we expect the 
description of gapped phases to be qualitatively correct over 
part of the phase diagram (although not the nature of the crit¬ 
ical points). In some parameter ranges, we find a good agree- 
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ment between the MF results and the more accurate DMRG 
calculations. 

To proceed we perform a MF decoupling of the interaction 
term H' in all channels respecting the symmetries of H. This 
leads to the MF Hamiltonian of the form 

= i 2 {TiOj/Bj + r2Pjaj+i + T\ajl3j+i + , 

(5.3) 

where r’s are the MF parameters. Note that, in the translation- 
ally invariant case, r' corresponds to the third neighbor hop¬ 
ping t' introduced in Eq. ( |3.3| ). These are determined based on 
the requirement that the ground state I'Fmf) of ^mf minimizes 
the expectation value {H)mv = ('I^mfI^I'I^mf)- This leads to a 
system of self-consistent MF equations for parameters r that 
read 


rameters (ti, T 2 , t^). This step must be performed numeri¬ 

cally. Second, we determine the topological phase character¬ 
izing the ground state of Hmv with these parameters. 

Because of the large parameter space involved in the anal¬ 
ysis it is instructive to start with the second step and enu¬ 
merate the possible topological phases of the noninteracting 
MF Hamiltonian ( |5.3| ). This is most easily done by studying 
its spectrum of excitations ( |531 ). We adopt a reasonable as¬ 
sumption that the gapped regions represent distinct topologi¬ 
cal phases with phase transitions marked by gap closings. We 
furthermore adopt t+ as our unit of energy and work with 3 
dimensionless coupling parameters 



T\ 


T 


h 




dEk ^ dEk\ 
5t2 


1 dEk 


(5.4) 

(5.5) 


and similar equations for T 2 and obtained by interchanging 
all indices 1 2. Here, N denotes the number of unit cells 

in the system, k extends over one half of the Brillouin zone 
(-TT, tt), which is now half the size of the BZ used in Sec. Ill 
because the unit cell has been doubled due to dimerization. 


Ek = 4^ a/IT+ sin 2 + y 


2 / h 

I T_ COS 2 ^ ^ 


(5.6) 

is the spectrum of excitations of Hmf and t+ = (ti + T2)/2. 

The physics of the MF approximation has a simple intuitive 
interpretation. In the noninteracting limit Hmf coincides with 
Hq. The interaction terms present in H' are seen to renor¬ 
malize the nearest neighbor hopping terms ti and t 2 via Eq. 


( |5.4[ ) and generate third-neighbor hoppings through Eq. ( |5.5| ). 
Second neighbor hoppings would violate T and are therefore 
not generated. For weak interactions it is easy to see that the 
MF equations imply an increase in hoppings with increasing 
|g| for g positive but decrease when g is negative. This obser¬ 
vation provides an intuitive explanation for the qualitatively 
different behavior of the model depending on the sign of g 
found in the previous Sections. We also note that the MF 
theory becomes exact in the strong coupling limit when, say 
g 2 = t 2 = 0 and for positive gi . [Notice that when g 2 = t 2 = 0, 
the Hamiltonian reduces to Eq. p.l0[ ), in which the operator 

= (^ 7 + 1)/2 is conserved. The ME approximation linearizes 
the fluctuations of hj and is exact if no such fluctuations are 
present.] In this limit Eqs. ( 5.4|5.5| ) reproduce the “ferromag¬ 
netic” ground state discussed in Sec. II.B. Eor negative gi the 
expected “antiferromagnetic” ground state breaks the symme¬ 
try under translation by two Majorana sites built into Hmf and 
is therefore not captured by this ME theory, although a more 
general ME theory could be constructed to describe this state. 

The ME analysis proceeds in two steps. Eirst, for given pa¬ 
rameters (/^i, ^ 2 ,,^ 1 ,^ 2 ) we And the ME Hamiltonian Hmf that 
best approximates H by solving Eqs. ( |5.4|5.5 ) to And ME pa- 


Eor Ek to be gapless both brackets under the square root in 
Eq.lj^ must separately vanish for the same momentum k. 
This imposes two conditions on three parameters (r, 5 -, s') and 
momentum k, leading to the conclusion that phase transitions 
occur at a set of two-dimensional surfaces in the 3-parameter 
space of Hmf- We expect this result to remain valid beyond the 
ME theory. When the translation symmetry is broken, a mass 
term becomes allowed in the low-energy effective Hamilto¬ 
nian ( |2.1| ). The mass m, then, is a function of 3 dimen¬ 
sionless parameters that can be constructed from couplings 
Phase transitions between massive phases cor¬ 
respond to m = 0 which imposes a single condition on 3 di¬ 
mensionless parameters, leading to the same conclusion as 
above. In the Luttinger Liquid -r Ising phase a mass term, 
imyRYL in the Ising sector and a pairing term A(ikRikL + h.c.) 
are allowed by symmetry in the dimerized model. A single 
condition is enough to make either m or A vanish correspond¬ 
ing to either a massless Ising model or a massless Luttinger 
Liquid. 

The first bracket in the spectrum in Eq. ( |5.6| ) vanishes (i) 
for k = 0 and all values of r or, (ii) for k = ko with r given by 

r=-^2-. (5.8) 

4sm^ (kol2)-3 


This solution exists only when r<-|orr>l. Now we must 
And for which values of 5- and 5’' the second bracket vanishes 
at these values of k. 

Consider flrst re (-|, 1). The second bracket vanishes for 
k = 0 when 5 '' = -5'. There is, therefore, a single phase tran¬ 
sition in the 5'-5'' plane indicated in Eig. [TT^. The transition 
takes place between a topological phase with v = 1 and a triv¬ 
ial v = 0 phase. This can be deduced by considering the limit 
r = 5-' = 0 in which Hmf coincides with the Hamiltonian of 
the Kitaev chain with nn hopping whose classiflcation is well 
known. Adiabatic continuity then insures that the invariant v 
remains unchanged unless we cross a phase boundary, hence 
the identiflcation of phases in Eig.pT^. 

Now consider r < - ^. In addition to the k = 0 solution we 
now have a solution at k = ko with ko given by Eq. {5^). The 
flrst solution implies a phase transition at s' = -s sls before 
while the second implies another phase transition line given 
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FIG. 11: Phase diagram of the mean-field Hamiltonian Hmf- Topo¬ 
logical phases for a) r G (-|, 1) and b) r < -1. The dashed phase 
boundary coincides with the solid line for r = -1, then rotates clock¬ 
wise with increasing |r|, eventually reaching slope -^1 as r ^ -oo. 
c) A chain with only third nn hoppings is equivalent to three decou¬ 
pled Kitaev chains. When t[ < the top and the bottom Kitaev 
chain will be in the topological phase while the middle chain will be 
trivial, leading to the invariant v = 2. In the opposite case the roles 
switch and we obtain y = -1. 


by 



The phase diagram is shown in Fig. \TT\ ). In addition to the 
y = 0,1 phases present before two new topological phases 
with y = -1,2 appear. The identification of these phases 
follows from observing that the gap closing described by Eq. 

) (dashed line in Fig. \n\ )) involves two Majorana modes 
and one thus expects y to change by +2 across this line. Al¬ 
ternately, one can consider the limit t+ ^ 0 (corresponding 
|r|,|v'| » \s\) in which the third nn hoppings dominate. In 
this limit the system breaks up into 3 weakly coupled Kitaev 
chains as illustrated in Fig. [HI- Two of these chains have an 
a operator on their left edge and one has a p operator there. 
The two possible phases in this limit are thus characterized by 
y = 2, -1. To verify this phase assignment we have also ex¬ 
plicitly computed the index y in a system with periodic bound¬ 
ary conditions and found agreement with Figs.[TT^,b. Similar 
analysis applies to the parameter region r > 1 with the same 
4 distinct phases possible but we find that this regime is never 
reached in the solution of MF equations ( |5.4|5.5| ) and is there¬ 
fore not relevant to our original interacting problem. 

We now proceed to analyze the MF equations ( |5.4|5.5| ). For 
concreteness and simplicity we set gi = g 2 = g, although it 
is no more difficult to analyze the general case. For g > 0 
we find that only two phases indicated in Fig. [TT^ appear in 
the MF theory and the transition occurs at ti = t 2 , as could be 
expected on the basis of symmetry when gi = g 2 . The tran¬ 


(5.9 



FIG. 12: Mean field parameters as a function of t 2 for fixed ti - gi - 
g 2 = 1- A jump at t 2 = 1 indicates the first order transition between 
topological phases with y = 0 and y = 1. 


sition is second order when g = 0 but becomes weakly first 
order for any g ^ 0. We have seen in Sec. II that the actual 
transition in the interacting problem remains second order up 
to large values of g so this is an artifact of the MF approxima¬ 
tion. We note, however, that for weak coupling the excitation 
gap at the transition point is exponentially small (~ so 
the MF theory provides at least a qualitatively correct descrip¬ 
tion of the transition at weak coupling (a spectrum with an 
exponentially small gap will be, for practical purposes, indis¬ 
tinguishable from a truly gapless spectrum). Fig. [T^ shows 
the discontinuity in the MF solution for t\ and T 2 for an in¬ 
termediate coupling strength g = 1, clearly indicating the first 
order transition. For larger g the discontinuity becomes more 
prominent but no new phases appear. 

For repulsive interactions g < 0 the behavior is qualita¬ 
tively different. Consider first the high-symmetry situation 
h = t 2 = t. Numerical solution of the MF equations ( 5.4|5.5| ) 
indicates that for g < 0 the effective nn hopping r decreases 
with increasing |g| from its initial value r = r at g = 0. At the 
same time r' also decreases, starting from zero and moving 
towards the negative values. When the velocity v = 2(t -h 3t') 
vanishes, a Fifshitz transition occurs with additional gapless 
branches of excitations making their appearance at ^ = ^o- To 
make contact with our work in Sec. Ill we plot in Fig. 4 the 
velocity v extracted from the MF solution as a function of t 
for g = -1. We observe that it vanishes at C - 3.41, rea¬ 
sonably close to the value 3.52 found from DMRG for a fully 
interacting problem. 

Now consider the gapped phases reached by perturbing the 
system away from the ti = t 2 symmetry line for g < 0. When 
t > tc such perturbations only give rise to y = 0,1 phases 
indicated in the phase diagram Fig. [HI- This is suggested 
by the solution of the MF equations and we have also veri¬ 
fied this by an exact numerical diagonalization of the full in¬ 
teracting Hamiltonian ( |5.1|5.2[ ). For t < MF theory sug¬ 
gests that all four distinct phases with y = -1,0,1,2 shown in 
Fig.[n]3 can be reached. This is in accord with the intuition 
that multiple species of low-energy gapless Majorana modes, 
once gapped by symmetry breaking perturbations, should give 
rise to phases with multiple MZMs bound to the edges of the 
chain. On the other hand, as indicated in Fig. [^c) we ex¬ 
pect the new y = -1,2 phases to be reached only when the 
effective third nn hoppings dominate. In the model defined by 
Eqs. ( |5.1|5.2 ) this will occur for relatively strong interaction 
strength because the third nn hoppings are absent when g = 0. 
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In this regime we expect the MF theory to be at best quali¬ 
tatively correct and our limited search using the exact diago- 
nalization method did not find conclusive evidence for these 
phases. The v = -1,2 phases can be stabilized by adding the 
third nn hoppings to the interacting Hamiltonian. They are 
allowed by symmetries and will generically be present in the 
physical system. We leave the detailed investigation of the 
resulting phases for future studies. 

We close this Section by noting that in the parameter re¬ 
gions where only v = 0,1 phases are present spontaneous 
dimerization in the geometry with open boundary conditions 
may favor the phase without MZMs. This can be seen by 
thinking about the strong coupling ground states of the model 
visualized in Figs. andClearly, the state with all MZMs 
combined into Dirac fermions will be lower in energy than the 
state with two unpaired MZMs at the edges. If there is a re¬ 
gion in the parameter space where only the v = -1,2 phases 
exist, presumably, by the same argument the v = -1 phase 
will be energetically favored over the v = 2 phase. 


VI. EXPERIMENTAL REALIZATIONS AND SIGNATURES 

As argued previously in Refs. (25] [26l all the ingredients 
are now in place to start experimentally exploring systems 
of strongly interacting Majoranas including the simple ID 
system discussed in this work. The most promising physi¬ 
cal platform is the superconducting surface of an STI where 
MZMs are bound in the cores of Abrikosov or Josephson vor¬ 
tices da- Experimentally, superconducting order has been 
induced in such surfaces by multiple groups and in several dif¬ 
ferent STI materials WfWSSj . The ability to tune the chemical 
potential to the vicinity of the Dirac point, required to bring in 
the regime of strong interactions, has also been demonstrated 
EIHSa. Recently, individual vortices have been imaged in 
these systems (541 and spectroscopic evidence indicative of 
MZMs in the cores of vortices has been reported (56l . ID 
structures such as those envisioned in this work most natu¬ 
rally arise in Josephson junctions built on the STI surface. 
When magnetic field is applied perpendicular to the surface of 
the STI a line of Josephson vortices is known to form inside 
the junction. Such Josephson vortices carry MZMs and their 
spacing can be conveniently controlled by the magnetic field 
amplitude. Evidence suggestive of MZMs in such devices has 
recently been reported i57l . 

The chain of ID interacting Majoranas could be realized 
in other physical systems that are known to host MZMs oi- 
El. This includes semiconductor quantum wires with strong 
spin-orbit coupling and the edge of a 2D topological insula¬ 
tor. If these are coupled to a periodic structure made of alter¬ 
nating superconducting and magnetic regions then MZMs are 
expected to form at the boundaries between the corresponding 
domains. In the presence of interactions the physics of such 
MZMs will be described by Hamiltonian ( |1.4| ) studied in this 
work. 

Scanning tunneling microscopy (STM) can provide valu¬ 
able experimental signatures of the phases and phase tran¬ 
sitions of this system when the model is realized by MZMs 


bound to vortex cores or other structures as discussed above. 
Importantly, the experimental ability to tunnel into a Joseph¬ 
son junction region has previously been demonstrated (58] 
|59l, both with STM and with planar contact tunneling. Be¬ 
low we address the characteristic signatures of various phases 
and phase transitions present in our model that are observable 
through the tunneling conductance. 

The tunneling current between the sample and the normal 
tip goes as (/) oc GR(-eV), where the retarded Green’s func¬ 
tion is (60l 

^oo 

GR(co) = -i\ (6.1) 

Jo 

Here the operator i/^o annihilates an electron at the tip. Due to a 
factorization of the time-ordered Green’s function into tip and 
sample correlators, the tunneling current captures the behavior 
of the temporal correlation functions of the Majorana chain, 
which are closely related to the spatial (equal-time) correlators 
at low energies due to the emergent Lorentz invariance. 

Several experimental signatures follow from this. For ex¬ 
ample, in the Ising phase, the Green’s function decays with 
the characteristic free-fermion exponent as x~^ and we find 
fising ^ V- At the tricritical point, the leading (with the small¬ 
est scaling dimension) fermionic operator has scaling dimen¬ 
sion 7/10. Therefore, as shown in Ref. (43l, the Green’s func¬ 
tion decays as and the tunneling current goes as 

/tci sgn(y)|yp/^ (6.2) 

In the 2-fold degenerate gapped (massive) phase, the power- 
law dependence of the tunneling current on V, which occurs in 
massless Ising and tricritical phases, becomes an exponential 
dependence. 

On the negative-g side, the Lifshitz transition is character¬ 
ized by a dynamical exponent z = 3, which changes the con¬ 
stant density of states of the Ising phase to a density of states 
at energy 6 proportional to 6“^/^. This leads to 

/Lifshitz ^ (6.3) 

Interestingly, at the Lifshitz transition the tunneling conduc¬ 
tance ^ diverges for small V. In the Ising-rLL phase, the 
Green’s function decays as x~^ to leading order, with sublead¬ 
ing corrections of the form gq p 21]) and 

Fig.iTjb)], indicating a tunneling current linear in V with sub¬ 
leading corrections (for small V) scaling as y(^+i/^)/2 

Although the contribution of the LL sector to the tunnelling 
current in the Ising-rLL phase is subdominant in V, it may 
be possible to see it clearly by doing STM near the end of a 
Majorana chain. There we expect a Luttinger-liquid contribu¬ 
tion that oscillates spatially at wave-vector 2ko while decay¬ 
ing with a ^-dependent power law with the distance from the 
end of the chain ED. Observing the values of ko = nj A and 
^ = 1/4, followed by a breakdown of the power-law depen¬ 
dence of the tunneling current on V, when entering the 4-fold 
degenerate gapped phase, can then signal the C-IC transition. 
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VII. CONCLUSIONS 

The spinful and spinless variants of the Hubbard model 
in ID, which respectively have on-site and nearest-neighbor 
interactions, are two widely studied canonical models of in¬ 
teracting Dirac fermions. Here, we studied a third canon¬ 
ical model, namely, the minimal ID model of interacting 
Majorana zero modes. These Majorana degrees of freedom 
have Hermitian creation operators and are effectively half of a 
(complex) Dirac fermion. In this case, the minimal interaction 
involves four sites. 

In light of the focused experimental efforts on realizing 
and manipulating Majorana zero modes as emergent particles 
in solid-state devices, the behavior of interacting many-body 
Majorana systems is of great theoretical as well as experimen¬ 
tal interest. We argued that ingredients necessary to begin 
experimental investigations of such systems with strong in¬ 
teractions are presently in place. Such investigations could 
naturally start with model defined in Eq. ( |1.4| ), this being the 
simplest interacting ID Hamiltonian one can construct with 
these degrees of freedom. 

Using a combination of analytical techniques based on 
field-theory and RG, mean-field calculations, and numeri¬ 
cal DMRG studies, we determined the full phase diagram 
of this novel strongly correlated system. The physics that 
emerges from this simplest model of interacting Majoranas 
is extremely rich and complex, revealing novel phases and 
phase transitions that are not present in the well studied Dirac 
counterparts. As previously shown in Ref. 14^ . this model 
provides one of the few examples of emergent space-time su¬ 
persymmetry for attractive interactions. In the present paper 
we extended the analysis of the model to the case of repulsive 
interactions. We found a novel z = 3 quantum critical point, at 
which a Lifshitz transition occurs changing the topology of the 
Fermi surface. On the weak-coupling (strong-coupling) side 
of the Lifshitz transition, we have one (three) species of low- 
energy Majoranas. In the Ising phase, the Majoranas are free, 
while on the strong-coupling side (the Ising-rLL phase), only 
one species remain free. The other two species combine into 
interacting Dirac fermions forming a Luttinger liquid, with an 
emergent Fermi momentum and particle number (despite the 
fact that particle number conservation is not a symmetry of 
the Hamiltonian). 

At even larger repulsive interactions, the Ising-rLL phase 
undergoes a transition to a 4-fold degenerate gapped phase. 
The nature of this novel phase transition, which fully gaps 
out a c = 3/2 CFT is reminiscent of the commensurate- 
incommensurate transition. As the emergent Fermi momen¬ 
tum varies with interactions, a term in the Hamiltonian, which 
generically has fast oscillations, becomes nonoscillatory at a 
particular commensurate wave vector and drives the transi¬ 
tion. 

When the translation symmetry of the chain is explicitly 
broken by dimerization, as often happens in ID systems, addi¬ 
tional phases can appear. We performed a brief survey of these 
dimerized phases using mean-field theory and found four dis¬ 
tinct gapped phases characterized by an integer topological 
invariant v. The latter takes values - 1,0,1,2, and can be inter¬ 


preted as the number of unpaired Majorana zero modes bound 
to the edge of the chain in the geometry with open boundary 
conditions. 

Our work extends the space of canonical ID models of in¬ 
teracting fermions from the Hubbard chain and its spinless 
variant to a third simple and experimentally relevant model 
describing the most natural interacting system composed of 
Majorana zero modes, revealing a plethora of novel phases 
and phase transitions. 

When this paper was almost complete, we became aware 
of (621, which contains related results. 
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Appendix A: Weak coupling 


While it is possible to diagonalize the noninteracting limit 
of Eq. ( |I.6| ) by a Bogliubov-DeCennes transformation, it is 
much simpler to Fourier transform the Majorana operators. 
For a chain of L Majorana operators, 7i, 72 , • • • 

yj = {-n<k< n) (Al) 


where k = ImIL with periodic boundary conditions or k = 
2n{n I/2)/L with anti-periodic boundary conditions. Invert¬ 
ing, we obtain: 


y{k) = 



This implies: 

{y{k),y{k’)) = dk-t, y(-k) = yHk). 


(A2) 


(A3) 


The Fourier transformed Hamiltonian is simply: 

H = 2t ^ y(-k)y(k) sin k. (A4) 

k 


This is already diagonalized and we see that we should iden¬ 
tify y{k) as an annihilation operator for k > 0 and as a creation 
operator for k < 0: 


y\k)y{k) sin^ - 

0<k<7T 


2 t 

sin(7r/L) 


(A5) 


In order to treat interactions using the renormalization group 
we focus on the low energy excitations, occurring at k near 
zero and tt where the dispersion relation becomes linear with 
slope V = 4t. In Schroedinger representation, we may write: 

yjit) K 2yR{vt - j) + i-iyiyiivl + j) (A6) 
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where 

jRivt-j) ^ + (A7) 

^ ^ 0<k<A 

rdvt + j) = ^ [e-‘'^^^‘^JM^-k) + e‘'^^^‘^jy{n-k)]. 

0<k<A 

Here A is a momentum cut-off, A 1. (We work in units 
where the lattice spacing is set to 1.) 7 /?/l(v^+ j) is a relativistic 
right/left moving Majorana fermion field. 


Appendix B: Spectrum of the Ising+LL phase 

Let us first review the finite-size spectrum for noninteract¬ 
ing relativistic Dirac fermions with Fermi wave-vector kp, 
with APBC. This depends on a number /, with |/| < 1/2, 
which is the fractional part of kpLIiln). The energy to add 
Nr fermions in the lowest energy states above the Fermi en¬ 
ergy on the right branch is: 

27rv 27rv /1 9 \ 

IT IT 

Taking Nr to be a negative integer, this formula also gives the 
energy to create Nr holes in the lowest energy states. Thus 
Eq. (|BTJ gives the energy for the lowest excited state with 
charge Nr relative to the charge of the ground state. After 
adding Nr particles in the lowest energy states we can make 
arbitrary particle-hole excitations. If NnR particles are raised 
by n energy levels the energy cost is (Inv / L)nNnR. So the 
energy for a general particle-hole excitation of right movers is 



FIG. 13: (a) The expectation value of pj, for PBC (computed with 
DMRG and retaining 700 states), goes to a nonzero value as ? ^ 0^ 
for g = 1 even in finite systems. As (pj) changes sign when t 
-t, this indicates a discontinuous jump in (pj). (b) Extrapolating 
{pj)\t^o+ to the thermodynamic limit strongly suggests that the jump 
survives at ^ ^ 00 . 


In addition, we may create excitations in the Ising sector. 
Adding Ir particles to the lowest energy states at small posi¬ 
tive k and Ip particles to the lowest energy states at k slightly 
larger than -n costs energy 


^El 


~T~ 


1 9 1 9 


(B6) 


Again 


- tX) 

27rv 

Er^ph = -j- ^NnRki. (B2) 

n=l 

The same formulas hold for left-moving excitations. Combin¬ 
ing them, it is convenient to define: 


N = Nr-\-Nl 

M = Nr-Nl. (B3) 


Ni = Ir + Il 

Ml = Ir-Il (B7) 


and Ni and M/ have the same parity. The formula for particle- 
hole excitations is more complicated in this case although it is 
well known. We won’t need it to analyze our DMRG data. It 
is convenient to shift N in Eq. ( |B5l > by twice the integer part 
of koL/ln, yielding Eq. (|3.10). 


Note that N and M must have the same parity: 

N = M, (mod 2). (B4) 


The excitation energy for an arbitrary low energy state is: 
2nv 


AE = 


1 

-(N - 2ff + -M^ + + N„L)n 

n=l 


(B5) 


These formulas carry over directly to the low energy exci¬ 
tations of the noninteracting Majorana system for t' < 
with kp replaced by k^. The analog of creating a right-moving 
hole corresponds to creating an excitation at k slightly less 
than -ko. 


Appendix C: Numerical evidence for the first-order transition 
between symmetry-broken phases for g > 0 


As discussed in Sec. II B, the first order transition corre¬ 
sponds to a splitting of the 4-fold degenerate ground states of 
Eig. 2 into 2-fold degenerate ground states plus a degenerate 
pair of excited states as t is turned on for g = 1. Depend¬ 
ing on the sign of t either the states with filled {{pj) = 1) or 
empty {{pj) = -1) Dirac levels are favored (first and third or 
second and fourth states in Eig. 2). With PBC, the surviving 
2-fold degenerate states have opposite fermion parity, as ar¬ 


gued above. [See Eq. ( |3.23| ] If we focus on the states with 
even fermion parity, 1st and 3rd in Eig. 2, then they are dis¬ 
tinguished by {pj). Eor r > 0 (^ < 0) the states with {pj) < 0 













16 


{{pj) > 0) occur. The transformation 727 ^ ~y 2 j^ which leads 
io t ^ -t, implies that {pj}(t) = -{pj}(-t). The first-order 
transition should be signaled by a jump in {pj} at t = 0. We 
found that this jump occurs even for finite systems. As seen in 


Fig. 13 a), the value of (p) saturates as t approaches 0^ (notice 
the logarithmic horizontal axis). In Fig.[^b), we extrapolate 


the value of {pj) for t ^ 0^ to £ 


• and find that it goes to 


0.70 + 0.01, clearly indicating a jump as we cross t = 0. 
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